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Abstract 



The magnetic and orbital orders for the bilayer manganites in the doping region 0.5 < 
X < 1.0 have been investigated from a model that incorporates the two orbit als at 
each Mn site, the inter-orbital Coulomb interaction and lattice distortions. The usual 
double exchange operates via the eg orbitals. It is shown that such a model reproduces 
much of the phase diagram recently obtained for the bilayer systems in this range of 
doping. The C-type phase with (vr, 0, vr) spin order seen by Ling et al. appears as a 
natural consequence of the layered geometry and is stabilised by the static distortions of 
the system. The orbital order is shown to drive the magnetic order while the anisotropic 
hopping across the eg orbitals, layered nature of the underlying structure and associated 
static distortions largely determine the orbital arrangements. 

PACS Nos. 75.47.Gk, 75.30.Et 

It has been realised [1, 2, 3] in the recent past that the physics of the region x > 0.5 is quite 
different from that in the x < 0.5 for the 3D manganites and one has to look at the heavily 
doped {x > 0.5) manganites from a different perspective. A similar situation prevails [4] in the 
bilayer manganites, the n = 2 member of the Ruddelsden-Popper series {R, A)n+iMnn03n+i 
(where R and A are rare-earth and alkaline-earth ions respectively) as well. The doping region 
< a; < 0.5 for bilayer manganites has been investigated in some detail and a rich variety 
of phases identified. These layered systems also show large magnetoresistance (MR) and a 
sequence of magnetic phases [5, 6] like their 3D counterparts. From a ferromagnetic (FM) state 
at low doping {x ~ 0) to canted antiferromagnetic (AFM) metallic to AFM insulating state 
between x = 0.37 to x = 0.48 have been reported [6, 7, 8]. At x = 0.5 there is a possible 
coexistence between charge ordered (CE-type) and layered A- type spin ordered state [7, 9]. 
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The region x > 0.5 has now been investigated [4, 10] carefully using neutron scattering and 
a succession of magnetic phases A— > C— > G has been observed. Between the A- and C-type 
phases (between 0.66 < x < 0.74), there appears a region of no well-defined long range order 
(LRO). Beyond x > 0.74 the AF C-type spin order is seen (along with a polytype, where the 
long c-axis is doubled). Interestingly, in the C-phase (or its polytype), the spins are aligned 
in the long basal plane b-axis, along which there is a distortion concomitant ai x — 0.74. In 
addition, both A-type and C-type phases have been found to be orbitally ordered. There is no 
evidence of canting of spins in the region x > 0.5. 

The role of orbitals on the underlying magnetic order is stressed [2] already in the context of 
the various magnetic structures of the 3D manganites. Models have been proposed [1, 3] for the 
manganites that incorporate the Cg orbitals and the anisotropic hopping between them. It was 
also realized that the inter-orbital interaction is quite crucial for the underlying orbital order 
[3, 11]. The use of such models to the bilayer manganites (like La2-2a;Sri+2xMn207 ) has only 
had limited success though [13, 14]. The quasi two-dimensional nature of the underlying lattice 
stabihses the A-type layered magnetic structure and the models have not been able to reproduce 
the observed C-type one-dimensional magnetic structure. The A-type AFM instability is indeed 
quite strong in the layered system (see fig. 1 in Ling et al [4]), extending from x = 0.42 to 
0.66. Moreover, at low temperatures, the CE-type spin and charge order seems to be absent and 
replaced by the A-type spin order [12], even at x = 0.5. On the other hand, there is a tetragonal 
to orthorhombic transition (elongation of the basal plane b-axis [4, 10]) near x — 0.74 where 
the C-phase appears. There is no buckhng of the octahedra associated with these distortions. 
The nature of spin and orbital ordering, as suggested by Ling et al. [4] and Qiu et al. [10], 
clearly points to the role of the electron-lattice coupling and the resulting elongation of the 
b-axis on the magnetic and orbital structure. Both the A- and C-phases are orbitally ordered 
and there is intimate connection between the preferred orbital orders, the lattice distortions 
and the magnetic order. 

The experimental observations and theoretical understanding generated for the heavily hole- 
doped 3D manganites quite naturally lead to a model for the bilayer manganites in the region of 
doping a; > 0.5 . The model incorporates the degenerate manifold and the physics of double 
exchange (DE) along with electron-electron and electron-lattice interactions. Such a model is 
given by 

H ^ JafYI Si-Sj -JhYI Si-Si - J2 4a,aCj,/3,(T + Hint + He-l (l) 

<ij> i <ij>cr,a,f3 

Here Si and Si represent the t2g and Cg spins at site i and Jh and Jaf are the Hund and super- 
exchange (SE) coupling respectively. The usual charge and spin dynamics of the conventional 
DE model operate here too, with additional degrees of freedom coming from the degenerate 
orbitals (a,/? take values 1 and 2 for the two e^, orbitals). The hopping across them is 
determined by the symmetry of Cg orbitals. The term Hint = U' 'fiiiani2a' describes on-site 
inter-orbital interaction. The intra-orbital term does not play a significant role for the typical 
values of Jh one is working with [11, 16]. The inter- bilayer exchange interaction is known to be 
at least a 100 times weaker [17] than the intra- bilayer one. Two bilayers are also well-separated 
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in an unit cell and intervened by the rare-earth ions. This allows us to consider only one bilayer 
for the calculations that follow. 

At X = 1 the Cg band has no electrons and the physics is governed entirely by the AF 
superexchange between the neighbouring t2g spins. On doping, the band begins to fill up (with 
nominal electron-density ^^). In the absence of electron- lattice coupling, the kinetic energy 
(KE) of electrons in the Cg band along with the attendant Hund's coupling between t2g and Cg 
spins begin to compete with the antiferromagnetic SE interaction leading to a rich variety of 
magnetic and orbital structures. The JT distortions, through the local electrostatic coupling 
(acting as an 'orbital magnetic field'), lift the degeneracy of the e^orbitals and affect the DE 
mechanism considerably. 

The coupling between the e^, manifold and lattice is incorporated through a term in H [18], 

He— I 9 ^ J '^i,mQi,m 
i,m 

where Qi^m{'^ = 1,2) are the even-parity local distortions of an MnOe octahedron and ri and 
T2 are the first and third Pauli matrices. A positive sign of g renders the 3^^ — orbital stable 
over — xp' orbital for distortion as there is negative charge on the surrounding oxygen 
ions. 

Writing Q^ x = ViSinOi and Qi^2 = TiCosOi, H^-i is diagonalised by the unitary transformation 

(cos— sin— \ 
. \ ^. . The choice of 9i 
S%TX'~^ COS j 

determines the orthogonal combination of orbitals and is dictated by the physics at hand. 
In addition, the orbital pseudospin operator turns out to be < 7^ >= (sin^j, 0, cos^j). The 
hopping matrices t^^fi along y, z directions, therefore, transform as Sit^'^'^S^. The rotational 
symmetry in the orbital space implies H^nt remains invariant. 

The diagonalisation of the KE part of H leads to two bands. In the pure (uncanted) phases 
the bands in A- and C-phases become purely two- and one-dimensional. However, even in 
the presence of canting there is little dispersion along the AFM aligned directions - a plane 
in C-phase or a line in A-phase. Typical values of the interaction and band parameters for 
the bilayer systems are in the same range as in the 3D manganites. The Hund coupling and 
Coulomb correlations are the largest scale of energy [16, 11] in the problem. Treating the t2g 

spins classically, the SE contribution to the ground state energy becomes EgE— ^2 ° {'^cos9xy + 
cos9z) where O^y and 9^ are the angle between the near- neighbour (nn) t2g spins in the xy plane 
and z direction respectively. 

For an uncanted homogeneous spin configuration in the ground state, we choose Sj = 
Soexp(iq.ri) where the choice of q determines different spin arrangements for the t2g spins [3]. 
We begin our discussion by considering the model without the Coulomb interaction terms U'. 
The nn Coulomb interaction and its effects will be dealt with later. 

Using the semi-classical approximation for the t2g spins the Hamiltonian (1) reduces to an 

8x8 matrix. The distortions are assumed to be uniform (r^ = \JQ\ + Ql = r). In almost 
all the manganites, the JT energy scales (2|5fr|) are nearly in the same order as the band- 
width, about leV . A typical value of \gr \ is therefore about 0.5eV" at a; = 0.55 [15], where the 
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tetragonal distortion is largest, Mn-Mn distance along c-direction shortest. The value of gr 
gradually decreases with increasing x as the c-axis elongates and vanishes by a; ~ 0.9. Around 
X — 0.75 there is a tetragonal to orthorhombic transition, with slight elongation of the basal 
b-axis disappearing by about x = .92. It is argued [10] that due to possible delocalisation of 
Cg electrons, the self-consistent JT scale around x — 0.55 could be much less. On the other 
hand there is evidence of charge ordering close to this region [8, 9, 10], which would lead to 
incipient localisation of charges. Nevertheless, the scale of static JT distortion used here is the 
bare value corresponding to an MnOg octahedron. 

We use mean-field approximation [3, 11, 14] to treat the Hamiltonian. This is shown to 
work quite well for the ground state properties [11] in the 3D manganites. The mean-field 
Hamiltonian is diagonalised at each k-point on a momentum grid. The ground state energy 
is calculated for different magnetic structures. We consider four different magnetic structures 
relevant for the experimental phase diagram. These are (with q values in the parentheses) 
A-type (0,0,7r), C-type (tt, vr, 0)-we call as C-type the usual C-phase with FM chains along 
c-direction, C'-type (tt, 0, tt) and the 3D AFM G-type (tt, tt, tt). The third one is the same as a 
C-type, only that its FM ordering is along b-direction as reported by Ling et al. The magnetic 
structure with minimum ground state energy is determined for each set of parameters {x, Jh, 
Jaf) for the range of doping (0.5 < ,t < 1) for a given distortion. Fig. la shows the ground 
state energy (all energies are measured in terms of = ^ — Q.2beV) with doping 0.5 < x < 1.0 
for typical values of exchange interactions for \gr\ = and 2.0 (along the c-axis) and Fig. lb 
shows the same with a distortion along b-direction nearly half the magnitude. The energies for 
\gr\ = are offset by 0.2 in order for better viewing. 

On shortening the bond lengths along c-axis, the energy of the C-phase rises while energies 
of both A and C phases go down. A-phase with its planar FM magnetic and orbital order 
(discussed below) is clearly favoured over the C-phase with out-of-plane FM magnetic (orbital) 
order. The C'-phase, with FM spin order along b-direction, also gains from the contraction in 
c-direction. This is even more apparent in Fig. lb where an elongation in basal b-direction 
stabilises C'-phase further. As reported in previous work [13, 14, 19] A-phase instability is quite 
strong in the layered manganites owing to the 2D structure of the DOS. The static distortion 
along b-direction stabilises both A and C phases, while the gain in stabihty of C phase is 
larger than that of A primarily due to its ID magnetic (and orbital) order along b-direction. 

The phases A and C are both orbitally ordered. Shown in fig. 2, the A-phase has planar 
a;^ — y'^ order while the C phase has 3y^ — order. The orbital densities do not change over 
continuously, there is an abrupt change across the A-C transition between the two sets of 
orthogonal orbitals indicating a first order transition between them. A strong orbital order is 
also seen [11] in exact diagonahsation study. Although the staggered orbital order is favoured 
close to half-filling, the second order t'^/\gr\ process is inoperative at this low electron-doped 
region where orbitals are mostly unoccupied. 

A phase diagram is then obtained in the \gr \ — x plane for typical values of JhSq and JafSq. 
It is observed (Fig. 3a) that with increasing \gr\ along c-direction, the C state stabilises slightly. 
The GC boundary is hardly affected as there are few electrons there. The large x part of the 
phase diagram is similar to 3D manganites primarily due to the absence of any significant 
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energy scales other than SE energy at such low electron-densities and reproduces the 3D AFM 
G-phase. The effect of elongation of the b-axis is more prominent as discussed above. The 
C'phase stabilises considerably over the A phase due to the changes in the occupied Cg DOS 
with enhanced orbital ordering of 3y^ — r^. The effect of change of bond lengths and consequent 
enhancement in bare hopping may stabilise A-phase somewhat when the c-axis contracts. The 
elongation in b-direction can also reduce the hopping in that direction thereby reducing the 
stability of C. The stabilisation coming from the static JT effects are expected to be stronger 
than changes coming from enhanced hopping at the doping regions considered. With changes 
in bond length less than 10% [4, 10], and the density of electrons low, this effect may not be 
large. In addition, the spin exchanges also depend on bond length (higher order in i as J ~ f^). 
Such effects are neglected in the presentation here. 

The phase diagrams in JhSq — x (Fig. 4) and JafSq — x (Fig. 5) reflect similar physics. 
To compare the theoretical phase diagram with experiments, in Fig. 5a, we have included the 
actual distortions between 0.5 < x < 0.92 with \gr \ = 2.0 at a; = 0.55 going down as x increases 
(by X — 0.75 the lattice nearly relaxes in the c-direction) [4]. The distortion in b-direction is 
smaller and occurs between 0.75 < x < 0.92. The phase diagram resembles the experimental 
one, albeit without the region of no spin order between 0.66 < x < 0.74. The C'phase in Fig. 
4a is more stable than that seen in experiments, covering this region of x where no apparent 
LRO is seen. Although the model recovers the C'phase seen in experiments, rather than the 
large A-type region observed in previous work [13, 14, 19], it overestimates the stability of this 
phase even without any static distortion. Note that there is a ferromagnetic phase in fig. 5 at 
very low J^p where the DE mechanism dominates. 

Canting of the magnetic structures Sj is included via Sj = So{sin(f)i, 0, cos0j) with 0j taking 
all values between and tt. In the G-phase, at large Jh, there is a small canting in the xy- 
plane (~ 8", inset in Fig. 4b), while 9z does not cant. The physics is quite similar to the 
3D manganites [3] and the x ~ region of bilayer systems [14]. At large in the G-phase 
the KE gain of the electrons through DE, via the generation of an FM component of the 
underlying t2g spins, more than offsets the 'cost' of tilting t2g spins away from magnetically ideal 
AFM state. Tilting in the xy-plane leads, of course, to a larger gain in KE than canting in 
9z, which remains insignificant. At smaller in the G-phase and in the A- and C'-phases, 
this mechanism is energetically inconsequential and we do not find any canting which is also 
reflected in the discontinuous (1st. order) change in the orbital order across A-C transition. 

We include the inter-orbital interaction term in the mean-field. As in 3D manganites [3], this 
term immediately stabilises C'phase. The ID instability of C'state is more favourably affected 
by the inter-orbital interaction and preferential occupation of orbitals due to U'. In addition, 
the higher electron-density in the A-phase makes this phase vulnerable to Coulomb interactions 
compared to the Cor G phase at lower electron-density [11]. 

The entire phase diagram with its magnetic and orbital order owes its origin to the compe- 
tition between DE mechanism, SE interaction, electron-lattice coupling and electron-electron 
interaction. In the region x ~ 0.5, where the electron-density is larger, the DE interaction 
via the degenerate Cg orbitals dominate. In the reduced dimensionality of the layered struc- 
ture, the planar dy.2_y2 orbital order along with DE coupling forces the ab-plane into an FM 
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configuration. The absence of long range correlation along c-direction and loss of tunnelling 
across the planes (driven by orbital order) induce AFM ordering in that direction and result 
in an A-phase. With a contraction of MuOq octahedra in the c-dircction, this phase further 
stabilises. Without a coherent charge transport in the c-direction in bilayer systems, the C- 
phase with (tt, tt, 0) magnetic order is unfavourable in comparison to the A-phase as already 
observed [13, 14]. Towards the x — 1 end, where the e^, levels are empty, the SE interaction 
brings about a (tt, tt, tt) magnetic order as in the 3D case. The C'-phase, on the other hand, 
allows for coherent tunnelling in the b-direction, its ID orbital order stabilises on contraction 
of the c-axis and elongation in the b-direction. At a certain x, as the electron-density reduces, 
this state stabilises over A-phase. The static JT distortions present in the system stabihses it 
until the SE interaction takes over at extreme low electron-doping. In the 3D manganites, the 
orbital order drives the magnetic order [1, 2] in the heavily hole-doped region. In the bilayer 
systems also, it is the orbital order, driven by the DE mechanism, anisotropic hopping across 
Cg orbitals and lattice distortions that seems to induce different magnetic phases. 

The scenario borne out here is markedly similar to the experimental phase diagram and 
orbital order (fig. 13 in Ling et al.) in the bilayer manganites. This also agrees quite well 
with the observed phases in Qiu et al. The existence of a region with no long range magnetic 
order aroimd x ~ 0.70 is quite possibly a result of the competing groimd states with such 
close energies (fig. 1). The A to C'transition being 1st. order in nature here there is a phase 
separated region (possibly dispersed due to long range Coulomb interactions). It would be 
interesting to look for inhomogeneous magnetic structures [10] or short range ordered phases 
(which are dispersed) in that region using more direct imaging techniques. It is also likely that 
with longer range Coulomb interactions included, charge ordered regions may stabilise close to 
X = 0.5, seen in several experiments [9, 10] recently. 
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Figure captions 



Fig. 1 Ground state energy of different magnetic phases versus hole concentration x > 0.5 with 
and without lattice distortions. In (a) the distrotion is in the c-direction while in (b) it 
is in the b-direction. The \gr\ — lines have been vertically offset by +0.2 to separate 
them from the lower bunch. 

Fig. 2. Orbital densities in (a) A- and (b) C'-phase at different values of parameters. In (a) 
the filled and open symbols are for dx2_y2 and d^^2_j.2 orbitals. In (b), they represent, 
respectively, the d^y2_^2 and dx2-y2 orbitals. The A- and C-phases are stable only in part 
of the range of x (see text). Note the sum of two orbital densities is equal to (1 — a;)/4, 
the actual electron density. 

Fig. 3. Magnetic phase diagram in \gr\ — x plane. Note the gradual shrinking of the A-phase in 
the region x > 0.5 while the G-phase remains nearly unaffected. 

Fig. 4. (a) Magnetic phase diagram in doping (x) - Jh^o plane is shown in solid line for exper- 
imentally relevant values of \gr\. In (b) is shown the effect of U' on the phase diagram 
(at \gr\ = 0). The solid hne is for U' — S and the dotted line in (a) and (b) are for 
\gr\ — U' — 0. In the inset in (b) is shown the canting of spins (away from tt) in the 
G-phase as a function of Jh^o- 

Fig. 5. Magnetic phase diagram in doping (x) - JafSq plane, (a) and (b) correspond to similar 
situations as in Fig. 4 (a),(b). 
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